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The two-phase issue in the 0(n) non-hnear cr-model: a Monte Carlo study 

B. Alles*** A. Buonanno*" and G. Cella'^ 

'^Dipartimento di Fisica and INFN, Piazza Torricelli 2, 56100 Pisa, Italy 

We have performed a high statistics Monte Carlo simulation to investigate whether the two-dimensional 0(n) 
non-linear sigma models are asymptotically free or they show a Kosterlitz-Thouless-like phase transition. We 
have calculated the mass gap and the magnetic susceptibility in the 0(8) model with standard action and the 
0(3) model with Symanzik action. Our results for 0(8) support the asymptotic freedom scenario. 
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1. INTRODUCTION 

The 2-dimensional 0{n) non-linear tr-model is 
defined by the action 
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together with the condition (/)(a;)^ — 1 for all 
spacetime points a;. In this equation (3 is the in- 
verse of the bare coupling constant. 

Perturbation theory (PT) predicts that this 
model is asymptotically free for n > 3. In partic- 
ular the exponential correlation length ^ on the 
lattice must scale as 
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where are the corrections to universal scaling. 
Here A = l/(n — 2). is a non-perturbative 
constant which for the standard action equals m 
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We define the magnetic susceptibility x the 
two-point correlation function at zero momen- 
tum. It scales as 
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where again is a non-perturbative constant. 
From equations (2) and (4) we conclude that in 
PT the ratio 
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tends to C-^/C"^ as we approach the continuum 
limit, (3 ^ oo. The corrections to asymptotic 
scaling dk depend on {afc} and {6^}. 

In a series of papers another scenario has 
been put forward for the model defined in (1). 
Under reasonable hypothesis the authors prove 
that there is no mass gap and that this model 
must undergo a Kosterlitz-Thouless-like {KT) 
phase transition at finite beta, Pkt- This implies 
that the ratio 
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should be constant as one approaches (3kt from 
below. Here 77 is a critical exponent. For the 0(2) 
model this exponent is 77 = 1/4. In |^ the authors 
show that the 0(3) model with the standard ac- 
tion on the lattice and r/ = 1/4 gives a constant 
for Rkt while the data for RpT displays a clear 
drop. 

Here we will show a progress report from an ex- 
tensive simulation performed on the 0(8) model 
with standard action and the 0(3) model with the 
tree- level improved Symanzik action If the 
constancy of Rkt for the 0(3) model is a gen- 
uine physical effect, then also for the Symanzik 
action we should see such a behaviour. The full 
account of our results with better statistics and 
using more corrections to asymptotic scaling can 
be found in j5|. 

2. SIMULATIONS 

In our simulations we have used the Wolff al- 
gorithm for the updatings as well as improved 
estimators to measure the correlation length and 
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Figure 1. The ratio Rkt for the 0(3) model with 
Symanzik action. 



magnetic susceptibihty. We have performed sev- 
eral millions of measurements for both quantities 
and verified the absence of autocorrelations. 

To calculate the correlation length we have 
measured the second moment The ratio 

^^^VC is less than few parts per mille, so within 
our statistical errors, we can use the formulae (2- 
3). 

We have chosen large enough lattice sizes L to 
keep finite-size effects under control. The ratio 
L/^ is 7 — 10. We have checked that these effects 
are few parts per mille. 

The largest systematic error comes from the de- 
viation from universal scaling of our data. These 
corrections are known up to 4 loops for the stan- 
dard action and up to 3 loops for the Symanzik 
action 0. 

3. RESULTS 

3.1. The 0(3) model 

In figure 1 we show the results for Rkt in the 
0(3) model with Symanzik action. They have 
better statistics than those of reference In 
constrast with our data are not constant. 

We do not show here (see |^) the data for Rpt- 
Again it is not constant although it displays bet- 
ter scaling than for the standard action 0. The 
fits for and agree with the prediction (3) 
and large-n calculations within 15 — 20%. 
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Figure 2. The non-perturbative constant for the 
correlation length. The 2, 3 and 4-loop results 
correspond to black circles, squares and diamonds 
respectively. The 2 and 3-loop results in the en- 
ergy scheme are the up and down triangles. 



Assuming finite-size scaling, it has been shown 
that this model presents asymptotic scaling start- 
ing from ^ « 10^5 g 

3.2. The 0(8) model 

In figure 2 we show the ratio between the non- 
perturbative constant 0|^^ as computed from our 
Monte Carlo data and the prediction (3) which for 
the 0(8) model is = 0.10544. If PT is correct 
and asymptotic scaling holds, this ratio should 
be equal to 1 (up to ~0.1 per mille because we 
measured the second moment ^^^•'). 

We show the data as obtained from the 2, 3 and 
4-loop approximation in eq. (2). The data in the 
scheme of the energy Q at 2 and 3-loop are also 
shown (we have used the energy measurements 
of reference jl^). All data seem to converge to- 
wards the PT prediction. A careful analysis of 
the next coefficients in the 1/n expansion sug- 
gests that further corrections should have small 
effects. The 4-loop data in this figure agree with 
(3) within 0.5%. 

In figure 3 the data for Rkt are shown. The 
data are far from constant. We show the results 
for two values of 77: 77 = 0.25 is the upper set of 
data (black circles) and 77 — 0.22 is the lower set 
(white circles). 
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Figure 3. The ratio Rkt for the 0(8) model. 
Black and white circles represent r] = 0.25 and 
rj — 0.22 respectively. The lines are the PT pre- 
dictions for these ratios. 
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Figure 4. The ratio Rpt for the 0(8) model. The 
successive orders correspond to circles, squares, 
diamonds and triangles respectively. 



The solid and dashed lines are the PT predic- 
tions for Rkt assuming = 0.103 (this is the 
value obtained from a best fit performed on our 
data for x; it agrees with large-n estimates within 
1%) and the prediction (3) for C{. We see that 
not only the curves are not constant but also that 
PT explains well its non-constancy. 

In figure 4 we show the data for Rpx- The up- 
per set of data is the lowest order prediction. Fur- 
ther corrections stabilize the result. The result 
clearly converges to a constant. Physical scaling 
is well reached at the largest correlation lengths. 
The corrections to universal scaling converge sur- 
prisingly well. 

The solid horizontal line is the PT prediction 
for the constant by using — 0.103 and the 
prediction (3) for the correlation length. We con- 
clude that our data are in fair agreement with 
PT. 

In conclusion our data do not support either 
KT or PT for the 0(3) model but they show clear 
agreement with PT for the 0(8) model. In the 
KT scenario one should explain why PT works 
so well for large n and large ratios L/£^. 



This work has benefited from many stimulating 
conversations with Andrea Pelissetto. 
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